Analysis of (K , K + ) inclusive spectrum with semiclassical distorted wave model 
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The inclusive K + momentum spectrum in the 12 C(K~ , K + ) reaction is calculated by the semiclas- 
sical distorted wave (SCDW) model, including the transition to the H~ bound state. The calculated 
spectra with the strength of the H~-nucleus potential —50, —20, and +10 MeV are compared with 
the experimental data measured at KEK with p K - — 1.65 GeV/c. The shape of the spectrum is 
reproduced by the calculation. Though the inclusive spectrum changes systematically depending 
on the potential strength, it is not possible to obtain a constraint on the potential from the present 
data. The calculated spectrum is found to have strong emission- angle dependence. We also investi- 
gate the incident K~ momentum dependence of the spectrum to see the effect of the Fermi motion 
of the target nucleons which is explicitly treated in the SCDW method. 

PACS numbers: 21.80.+a, 25.80.Nv, 24.10.Eq, 24.50.+g 



I. INTRODUCTION 

Hypernuclear physics has been disclosing characteris- 
tic features of strong interactions in baryon systems with 
strangeness. The A-iV interaction has been fairly well 
studied through analyses of experimental data on from 
light to heavy A hypernuclei. E-JV and S-7V interactions, 
however, are not well explored because of insufficient ex- 
perimental data. In particular, the 3-iV interaction is 
much ambiguous for this reason despite the fact that its 
knowledge is the most important for understanding the 
physics of octet baryons as a probe of the strangeness 
S = — 2 sector. 

Since the information from direct Y-N scattering ex- 
periments is extremely limited, the data on the produc- 
tion of a hyperon Y in strangeness exchange reactions 
(tt + ,K + ), (it~,K + ), (K~,K + ) etc. have been a princi- 
pal source of information on the single-particle potential 
Uy felt by Y in the final hypernucleus that reflects the 
properties of the underlying Y-N interaction. The well 
depth of the 2-nucleus potential Us was conjectured to 
be about 14 MeV by a distorted wave impulse approxima- 
tion (DWIA) analysis of the inclusive K + spectra in the 
threshold region of the (K~ , K + ) reaction [lj,|2|. This 
potential depth is not conclusive, however, because of 
the limited statistics and resolution. The situation will 
be improved in the near future, since the (K~ , K + ) ex- 
periments with high-intensity K~ beams with a better 
resolution are prepared in the J-PARC project at KEK. 

In the recent papers 0, , the inclusive spectra of the 
(ir~ , K + ) X formation reaction were analysed by means 
of the semiclassical distorted wave (SCDW) model [5| and 
information on the S-nucleus potential was extracted. 
In this paper, we calculate, using SCDW, the inclusive 
K + spectrum in the (K~ , K + ) S~ production reaction 
in an attempt to determine the S~-nucleus potential. In 
view of the on-going experimental progress, we try to 
elucidate the angle dependence and the incident energy 
dependence of the inclusive cross section. The SCDW 
model has been successful in quantitatively describing 



(p, p'x) and (p, nx) inclusive cross sections for wide range 
of energy transfers and emission angles without any free 
adjustable parameter 0,00- The single free parame- 
ter in the present SCDW calculation of the strangeness 
exchange reactions is the strength of the S~-nucleus po- 
tential Uy- 

The formulation of the SCDW model for describing the 
(K~, K + ) spectrum is described in Sec. |nj for quasifree 
region in Sec. Ill Al and for bound state region in Sec. 
I11BI In Sec. 11111 descriptions are given of the kaon dis- 
torted waves, the single-particle states of target nucleons 
and the S~ hyperon, and the strength of the elemen- 
tary process K~ + p — > K + + E - . Numerical results of 
the calculations for the (K~,K + ) spectra are compared 
with the experimental data and discussed in Sec. IIVI A 
summary and conclusions are given in Sec. [3 



II. FORMALISM 

A. SCDW for inclusive spectrum 

The inclusive double differential cross section in the 
center of mass (cm.) system for the 12 C(lf~, K + ) reac- 
tion with production of a S - hyperon is given in the first 
order DWIA by 
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where -E c ,rcd = E c E c ,a/ {E c + E c ^a) is the reduced energy, 
p c is the relative momentum between the kaon and the 

target 12 C nucleus, Xc + '°'' is the distorted wave of the 
kaon at energy E c for c = i and / where i and / stand 
for the initial and the final state, respectively, and E^a 
(Ef,A) is the total energy of the target nucleus (hyper- 
nucleus) in the cm. system. The factor AEiEf(2n) 6 is 
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the product of the normalization factors of the \ c given 
by 



(Xc(p c )\Xc(p' c )) = 2E c (2ir) 3 5 3 (p c - p' c ) 



(2) 



The suffices a and f3 denote the states of the nucleon 
hole and the unobserved S _ hyperon in the final hyper- 
nucleus, respectively, and <j> a and <f>p are the correspond- 
ing single particle wave functions with energies e a and 
Ep. The delta function guarantees energy conservation 
where u> — Ei — Ef is the energy transfer. The i-matrix 
in coordinate representation of the elementary process 
K ~~ + p — » K + + S~ is denoted by v that depends on 
energies and momenta of the particles in the reaction. In 
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the present calculation we assume the range of v to be 
zero. 

In the laboratory system, in which most experimental 
data are given, the cross section is given by 
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where L denotes the laboratory system and J = 
d{E f ,Q, f )/d{p),n J j : ) is the Jacobian. 

Expanding the squared modulus in Eq. (JJJ, one ob- 
tains 
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where the nonlocal kernel K(r,r') is given by 
K(r,r') = Y,^ r ^(r)4>p{r'W a {r') 
xS(sp - e a -uj). 



(5) 



The summation over a extends to all the hole states 
below the Fermi level. The summation over (3 is unre- 
stricted since there is no Pauli principle for the produced 
hyperon. 

For the summation over a, we use the Wigner trans- 
form 
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where R = r+ J" and s = r' — r denote the midpoint 
and the relative coordinates of r and r', respectively. 
The summation over (3 is replaced by the integrals 
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over the outgoing direction Qp and the energy Ep of the 
produced hyperon in state (3 in the asymptotic region, 
and fiy is the reduced mass of the S _ -nucleus system in 
the final state. 

In the SCDW model, the local semiclassical approxi- 
mation (LSCA) is used for the distorted waves %c KJ 



Xc (R±k s ) = Xc(R)e ±ik ^ R y^ S , (8) 
for small s/2, where hk c (R) is the local momentum. The 



direction of k c (R) is calculated by 

Re[ x ( ±) *( J R)(-i)V X ( ±) (H)] 
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and the magnitude is determined by the energy- 
momentum relation at R, 

h 2 c 2 k 2 (R) + 2E c , Icd (U R {R) + Vcoui(-H)) 

-V 2 oul (R) = El - m 2 K c A , (10) 

where ttik is the kaon mass 494 MeV/c 2 , Vcoui(-R) is the 
Coulomb potential, and Ur(R) is the real part of the 
distorting potential for Xc with E c . 

The LSCA for the distorted waves \c is expected to 
work well, because the nonlocal kernel K(r,r'), which 
can now be written as 



K(R— -a, R + ~a) = 
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is appreciable only for small s because of the complete- 
ness Y,p<l>p{R - s/2)4>f}[R + s/2) = 5(s). It is noted 
that the range of the kernel is shorter than in the cases 
of (p,p'x) and (p, nx), in which the Pauli principle re- 
stricts the summation over (3. 

In the present calculation, we make use of the LSCA for 
the unbound wave functions <pp in the S _ -nucleus poten- 
tial Us, hence, (j)p(R±s/2) = <j>p(R)exp[±ikp(R) - s/2] 
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with the direction given from the flux of 4>p{R) as 
in Eq. |JS} and the magnitude fc /3 (J?) = {2\iy\sp — 

UsiR)]/^ 2 } 1 ^ 2 , where Us(R) is assumed to be real as 
described in Sec. II I II The approximation is valid when 
ep is large compared with Us(R) and kp{R}j is a slowly 



function of R. 

Using the approximations described above, the dou- 
ble differential cross section for the inclusive reaction be- 
comes 
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where the first delta function on the right hand side implies the local momentum conservation. 



(12) 



B. SCDW for exclusive spectrum 

We apply the SCDW model also to the (K~,K + ) spectrum in the region of bound states of hypernucleus. The 
differential cross section in the laboratory system for the transition to a S~- n B bound state is given by 
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where J' = dil / / d&j is the Jacobian. Compared with the inclusive spectrum, this formula has neither the energy 
delta function nor the summation of a or (3. Corresponding to a transition from proton-hole state a' to a S~-particle 
state (3 1 with the energies e a ' and sp>, respectively, there should be a spike at the point with the energy transfer 
u = ep> — e a i. 

We use the Wigner transforms of the product of wave functions of the nucleon hole and those of the hyperon: 

-ik„,s 



<Mr)C,(r') = J dk a - $ al {R,k a ,)e- lK «' s , (14) 
^(r)^(r') = J dkp,<$>p,{R,kp,y k ?' s . (15) 
Using Eqs. (|14|) and i|15|) and changing the variables of integration to (R, s), one can write Eq. (|13fl as 
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In this case, the nonlocal kernel K(r,r') = ^) a >(r)^,(r')4>pi(r)4>p>(r') has the range s of the order of the target 
nuclear radius. Since we still use the LSCA for the \c, the accuracy of Eq. I|16|l depends on that of the LSCA with 
such large s, which will be discussed in Sec. IIV Bl 



III. WAVE FUNCTIONS AND TRANSITION 
STRENGTH 



system reads 

[h 2 c 2 V 2 r + {p\fc 2 - 2E^(UHr) + U c L oul (r)) 



+ (V£ oul (r)) 2 ]x L c (r) = 0, (17) 

The distorted waves of the K ± are calculated as so- where r is the displacement of the kaon from the center 
lutions of the standard Klein-Gordon equation with dis- of mass of the target nucleus, p\ (E^) is the momen- 
torting potentials U c . The equation in the laboratory turn (energy) of the kaon in the laboratory system, and 
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U\ and Vboui are, respectively, the distorting and the 
Coulomb potentials defined in the laboratory system for 
the states c — i and /. In the present calculation, we ap- 
proximate Eq. 1171) by the equation in the cm. system 
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replacing p\, E\, £/ c L , and Vboui by p c , E c ^ cd , U c , and 
Vboui, respectively, in order to include the nuclear recoil 
effect. We use distorting potentials of the standard tp 
approximation of Ref. 8] : 
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where = M c ^aWkn/M(E c + E C: a) is a kinematical 
factor resulting from the transformation of the KN scat- 
tering amplitude from the K-N to the AT-nucleus cm. 
systems, <jkn is the isospin-averaged KN total cross sec- 
tion, /.R-jv(O) is the KN forward scattering amplitude in 
the free space, and p{r) is the nucleon density distribu- 
tion of the nucleus in the A"-nucleus cm. system. The 
R e [/xw(0)] and <Jkn are obtained as the isospin- average 
of those given in Ref. in parameterized forms. The 
masses of the nucleon, the target nucleus, and the hyper- 
nucleus are denoted by M, Mi y A, and Mf,A, respectively, 
and Wkn ikKN) denotes the total kaon- nucleon energy 
(wave number) in their cm. system. 

Actually, the U c do not agree with the optical poten- 
tials that fit the experimental total cross sections ctt 
of K±. Fi 

gure^shows t ne °t for K^ on 12 C calculated 
with the (dashed line) and the U c (solid line). It turns 
out, however, that the inclusive (K~, K + ) cross section 
calculated with the U c agrees with that calculated with 
the within 10% which is within the experimental er- 
ror of the (K~ , K + ) spectrum. We use, therefore, the 
potentials U c in our calculations. 

The single-particle wave functions and the density dis- 
tribution function p(r) of the target nucleus are calcu- 
lated in the Hartree-Fock model of Campi and Sprung 
|ll| . The unobserved hyperon S - is assumed to propa- 
gate in a Woods-Saxon potential, 
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1 + cxp{(r - R)/a} 

with R = 1.1A 1 / 3 fm and a — 0.65 fm. The calculation 
is performed with a series of Vs = —50, —20, and +10 
MeV. The flux of the S~ decreases in the propagation 
because of its interaction with other nucleons. We take 
that effect and also the experimental energy resolution 
into account by means of convolution of the calculated 
spectra with a Lorentz form of width T. 

The transition strength \v\ of the elementary process 
K~ + p — > K + + S~ is assumed to be related to the 
differential cross section in the free space 

da _ 1 E p E s k K + 2 
dn ~ (4tt) 2 W 2 fc^~ 
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1.5 2 
P K [GeV/c] 

FIG. 1: The total cross sections of the K~ and the K + colli- 
sions on 12 C as functions of the incident momenta of K~ and 
K + in the laboratory frame. The optical model values with 
the tp potentials U c and the optical potentials U are rep- 
resented by the solid lines and the dashed lines, respectively, 
and compared with the experimental data of Ref. 10j • 



in the on-shell approximation, where W is the invariant 
energy, and E p (E~,) is the energy of the proton (the 
S~ hyperon) and kx- (kx+) is the momentum of K~ 
(K + ) in the K~-p (K + -E.~) cm. system. We use the 
elementary cross section parameterized in Ref. | 12| as a 
function of s and the scattering angle 9. In our model 
described in the previous section, W and 9 are given by 
the local momenta of the particles K~ and p in the initial 
state and K + and 2~ in the final state at each collision 
point R. 



IV. RESULTS AND DISCUSSION 
A. Transition to continuum states 

In Fig. [21 the momentum spectrum of K + in 
12 C(A~, K + ) at incident momentum p K - = 1.65 GeV/c 
calculated with Vs — —50 (dotted line), —20 (solid line), 
and +10 MeV (dashed line) are compared with the KEK 
experimental data [l3|. Since the data are average over 
the emission angles 9 K + between 1.7° and 13.6°, the cal- 
culated values are correspondingly averaged. The spec- 
trum shown is after the Lorentzian convolution men- 
tioned in the preceding section. The result turns out 
to be quite insensitive to the width T of the convolution. 
In Fig. [5] the half width T/2 is tentatively taken to be 
2 MeV, the same value as for the transitions to bound 
states discussed later. 

The shape of the spectrum is reproduced by the cal- 
culation, but the absolute magnitude of the cross section 
is underestimated by about a factor of 2. A possible 
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2: Momentum spectrum of K + in 12 C(K ~ , K + ) at 
= 1.65 GeV/c averaged over the emission angles between 
= 1.7° and 13.6°. The data of Ref. [l3| are compared 
with theoretical values calculated with Vs = —50, —20, and 
+ 10 MeV represented by the curves as indicated. 
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FIG. 3: Emission angle dependence of the momentum spec- 
trum in 12 C{K~, K+) at p K - = 1.65 GeV/c with V s = -20 
MeV. The spectra at 6 K + — 0° and 13° are shown together 
with the average over 1.7° < 8 K + < 13.6°. 



reason on the theoretical side is that we do not correct 
the lack of the center-of-mass motion of the target wave 
function. Since the momentum transfer from K~ to K + 
is large, the effect is not negligible. This problem in the 
context of the SCDW method has yet to be clarified. On 
the physics side, we ignore multi-step processes. Their 
contributions have been estimated by Nara et al. [T^ 
by means of the intranuclear cascade model (INC), and 
shown to be small in the case of 12 C(K~, K + ). Since 
INC is a classical mechanical simulation, however, rein- 
vestigation by means of the quantum mechanical SCDW 
model is desirable, even though the magnitude of the 
cross sections in Fig. [5] is almost the same as that of the 
one-step process in given by the INC calculation of Ref. 
[l2|. Another possible reason for the discrepancy is the 
in-medium modification of the elementary process which 
is not taken into account in the calculation. There is 
also experimental uncertainty of about 30 % in the mag- 
nitude of the elementary cross section. Figure [5] shows 
that the V^-dependence of the angle-averaged spectrum 
is not very strong. The peak position shifts naturally 
as the potential strength changes. It is obvious, how- 
ever, that one hardly determines the value of Vs from 
the present experimental data. 

Figure[3]shows the same calculated momentum spectra 
as in Fig. Hat 8 K + = 0° (dotted line) and 13° (dashed 
line). One sees that the 9k+ dependence is very strong. 
This is in contrast to the almost same spectra at 0° and 
13.6° in Ref. [l4( obtained with angle-averaged elemen- 
tary cross sections. This clearly shows the significance of 
the use of angular dependent elementary cross sections 
in the present SCDW model calculation. 

Figure 01 shows the dependence of the calculated 
12 C(K ~, K + ) energy spectrum on the incident momen- 



tum pk-, which is not very strong. The dotted, solid, 
and dashed curves denote the calculated energy spectra 
at p K - — 1.50, 1.65, and 1.80 GeV/c, respectively. It is 
interesting to note, however, that the change of the po- 
sition and the hight of the peak with the change of Pk- 
from 1.50 to 1.80 GeV/c is not simple, unlike that of the 
elementary cross section of Ref. 01 used in the present 
calculation. Analysis of the calculation shows that the 
Fermi motion of the target nucleon is responsible for this 
somewhat complicated Pk- dependence of the spectrum. 



B. Transition to bound states 

If the potential Us is attractive and sufficiently strong, 
the hypernuclear system S - - n B can have bound states. 
Table|l]shows the energy eigenvalues of the single-particle 
states of S~ in Us for the depth of Vs — —20 and —50 
MeV. Transitions to states with a nucleon hole in the 
0s7, 2 -orbit are treated with the convolution of the calcu- 
lated spectra with the wide width of 9.2 MeV [l5| . DWIA 
calculations of the (K ~, K + ) cross sections are made as 
in the case of transitions to continuous states (TCS), ex- 
cept that the LSCA is not used to the S~ bound state 
wave functions. Figure [S] shows the calculated momen- 
tum spectra of K + at p K - = 1.65 GeV/c and 9 K + = 8° 
in the region of p^ + where the peaks corresponding to 
the bound states appear. The calculated spectrum is con- 
voluted by Lorentzian with the half width T/2 = 2 MeV 
which is the sum of 1 MeV due to the decay — ► AA 
and 1 MeV due to the resolution of the detector in the 
experiment. 

In Fig. [31 the dotted, solid, and dashed lines corre- 
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FIG. 4: The energy spectrum of K + at 6 K + = 0° at p K - — 
1.50, 1.65, and 1.80 GeV/c calculated with V~ = -20 MeV 
and r/2 = 0, are plotted as functions of Mhy — Ma where 
Mhy and Ma are the masses of the hypernucleus and the 
target nucleus, respectively. 



TABLE I: Eigenstates and eigenenergies of H in the single- 
particle potential Us- 



V s [MeV] 


E~ [MeV] 




Os-state Op-state ls-state 


-20 


-8.3 -0.6 


-50 


-28.7 -13.0 -2.6 



spond to the calculated spectra with Vs = —50, —20, 
and +10 MeV, respectively. One sees two (three) peaks 
corresponding to the transitions from the Op^-state to 
the 0s- and Op-states (0s-, Op- and ls-states) in Table H] 
for Vs = -20 (-50) MeV, while for Vs = +10 MeV, as 
anticipated, there is no peak. 

An essential difference from the case of TCS is that 
there is no summation over the final states and, conse- 
quently, no short ranged kernel discussed in Sec. ITTlthat 
would restrict the distance of the source points of the 
outgoing distorted waves \f that can interfere with each 
other. Therefore, one needs the LSCA extrapolation of 
Xf over a long distance, the validity of which needs to 
be examined. Numerical tests show that the accuracy 
of the LSCA extrapolation from a point on the nuclear 
surface to another point across the nuclear diameter is 
within about a factor of 2. This is the sort of accuracy 
of the calculations shown in Fig. [5] 

The presence of a bound state peak in the spectrum 
shows that the S~-nucleus interaction is attractive and 
strong enough to support that particular bound state. 
Through the model calculation, one can estimate the 
strength of the peak that can suggest the experimen- 
tal precision required for measuring it. In this regard, 
we consider what the results of the present work suggest 



FIG. 5: The 12 C(K~,K + ) spectra near the S"- n B hypernu- 
clear bound state region at p K ~ = 1.65GeV/c and 6 K + = 8° 
calculated with Vs = -50, -20, and +10 MeV. 



is significant despite the numerical ambiguities discussed 
above. 



V. SUMMARY AND CONCLUSIONS 

The inclusive K + momentum spectrum in the 
12 C(K ~, K + ) reaction is calculated by the semiclassical 
distorted wave (SCDW) model. The same method of the 
calculation is used as in our previous works on (7r _ , K + ) 
and (ir + 1 K + ) inclusive spectra (3,3- The transition to 
the S~ bound state is also considered. The accuracy 
of the local semiclassical approximation (LSCA) for the 
discrete part of the spectrum is not so good as for the 
continuous part, but the error is estimated to be within 
a factor of 2. The standard tp potentials are used for the 
if ± - 12 C distorting potentials. Although the total cross 
sections of the K~ and K + scatterings on 12 C calculated 
with those potentials do not reproduce the experimen- 
tal data very well, the resulting (K~ ,K + ) spectra agree 
with the ones calculated with the optical potentials U ± 
within about 10 %. In order to improve a model for kaon 
distorted waves, one needs more experimental data in the 
momentum region of interest. 

Since we can explicitly take into account the depen- 
dences of the elementary K~ + p — > K + + S~ cross 
section on the invariant mass squared and the scatter- 
ing angle, the present calculation is more realistic than 
those, e.g. in Ref. ^4[, in which some constant aver- 
age value of the elementary cross section is employed. 
The S~-nucleus potential is assumed to be of a Woods- 
Saxon form. The inclusive momentum spectra with the 
strength of the 12 C(K~ , K + ) reaction calculated with the 
strength of the potential Vs = -50, -20, and +10 MeV 
are compared with the experimental data at p K - = 1.65 
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GeV/c. The shape of the spectrum is reproduced by the 
calculation. However, the magnitude of the cross section 
is underestimated by about a factor of 2. Though the 
inclusive spectrum changes systematically depending on 
the strength of V3, it is not possible to obtain a con- 
straint on Vs from the present analysis on the basis of 
the available data. 

There are several possible reasons for the underestima- 
tion. There should be the contribution of the multi-step 
processes to (K~,K + ). We will investigate, in the next 
step, the influence of the multi-step processes by means 
of the SCDW model for these processes 0|. The effect of 
in-medium modifications of the elementary process may 
not be negligible, although it can reduce the cross sec- 
tion. Another source of the discrepancy is experimental 
uncertainty of about 30 % in the elementary cross sec- 
tions used in the calculations. It is obvious that we need 
experimental data of the (K~, K + ) spectrum and the el- 
ementary cross sections with better accuracy for a more 
quantitatively reliable analysis. 

The calculated spectrum is found to have strong 
emission- angle dependence. When analyzing (K~,K + ) 
experimental spectra with angle- averaging involved, it 



is important to be aware of such dependence. It is 
seen from the spectra calculated at incident momenta of 
Pk- = 1.50, 1.65, and 1.80 GeV/c that the cross section 
at the lower side of the excitation energy hardly depends 
on the K~ incident momentum. The close examination 
of the calculation shows that this consequence is due to 
the Fermi motion of the target nucleons, which is prop- 
erly treated in our SCDW framework. The discrete part 
of the spectrum of K + will be useful for extracting in- 
formation on the S~-nucleus potential from the future 
(K~,K + ) experimental data, because the possibility to 
detect a peak structure below the threshold has not been 
ruled out yet experimentally. 
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